A density functional theory for general hard-core lattice gases 
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We put forward a general procedure to obtain an approximate free energy density functional 
for any hard-core lattice gas, regardless of the shape of the particles, the underlying lattice or 
the dimension of the system. The procedure is conceptually very simple and recovers effortlessly 
previous results for some particular systems. Also, the obtained density functionals belong to the 
class of fundamental measure functionals and, therefore, are always consistent through dimensional 
reduction. We discuss possible extensions of this method to account for attractive lattice models. 
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Despite the crucial role that lattice models have had in 
the development of Statistical Physics, when one looks for 
such models in the literature of density functional theory, 
the results are scarce. In the last few years, though, some 
of the most classical approximations have been extended 
to lattice systems 0, H, 01 and used to study different 
phenomena (like freezing and fluid-solid interfaces 0, Q 
or confined fluids Q)- Also very recently, fundamental 
measure (FM) theory has been added to the list through 
its formulation for systems of parallel hard hypercubes in 
hypercubic lattices The construction mimics that of 
its continuum counterpart , and like the FM functional 
for hard spheres, it is obtained from a zero-dimensional 
(Od) functional (a functional for cavities holding one par- 
ticle at most). This theory possesses a remarkable prop- 
erty: dimensional crossover, which allows obtaining the 
functional for d — \ dimensions from the one for d di- 
mensions by confining the system through an external 
field to lie in a (d — l)-dimensional slit. Dimensional 
crossover has been applied to the already mentioned sys- 
tem of parallel hard hypercubes in order to obtain FM 
functionals for nearest-neighbor exclusion lattice gases 
in two-dimensional (square and triangular) and three- 
dimensional (simple, body-centered and face-centered cu- 
bic) lattices p. 

This increasing interest in density functionals for lat- 
tice models has several motivations. On the one hand, 
some systems are particularly difficult to study using con- 
tinuum models. For them, lattice models provide conve- 
nient simplifications. This is the case of glasses or 
fluids in porous media to name only two. On the 
other hand, lattice models cover a wider range of prob- 
lems, many of which do not even belong to the theory 
of fluids (like roughening or DNA denaturation [lo{. 
to name only two) and so have never been studied with 
density functional theory. Finally, from a purely theo- 
retical point of view, these extensions are also interest- 
ing because they reveal features of the structure of the 
approximate functionals which are hidden or at least not 
apparent in their continuum counterparts (this is the case 



of FM functionals). 

In this letter we propose a simple systematic procedure 
to construct a FM functional for any hard-core lattice 
model. The construction is based on the dimensional 
crossover of this theory, much like the latest versions for 
continuum models. 

Let us begin by realizing that all FM functionals for 
lattice models studied in Refs. 0, share a common 
pattern, namely, the excess free energy can be written as 

(1) 
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where C denotes the lattice, I is a set of indices suit- 
ably chosen to denote the different weighted densities 
n'^^\s) = X]tGCfc(s) /°(*)' are integer coefficients which 
depend on the specific model, $o(??) — V+{^~v) ln(l — ?y) 
is the excess free energy of a Od cavity with average occu- 
pancy < rj < 1, p{s) is the density profile of the system 
(specifically, the occupancy probability of node s) and 
Cfe(s) is, for each k € T, a finite labeled subgraph of the 
lattice placed at node s (vertices are labeled with node 
vectors). The shape of the graphs Cfe(s) also depends on 
the model. From the definition, n('^^(s) appears as the 
mean occupancy of the lattice region defined by Cfc(s). 

For the sake of clarity we will illustrate this formal 
setup and the arguments to come with a simple exam- 
ple: the two-dimensional square lattice gas with first and 
second neighbor exclusion. With the help of a diagram- 
matic notation already introduced in jl3l | , the excess free- 
energy functional for this model takes the form 

P^'^'ip] = E (O) - *0 («) - $0 il) + ^0 (o)] , 

(2) 

where the diagrams represent the (four in this case) 
weighted densities H = n'^^'^^s), = n(i^")(s), J = 
n(0'i)(s) and o = n(0'0)(s), where [s = (31,82)] 

n('=-^-^)(s) = EEpG^i + V'*2+j). (3) 

i=0 j=0 
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FIG. 1: Examples of Od profiles corresponding to maximal 
cavities for the first and second neighbor exclusion lattice gas 
in the square lattice (a) and the nearest neighbor exclusion 
lattice gas in the triangular lattice (b). The density can only 
have nonzero values at the black nodes (which define maximal 
cavities)). 

This notation uses explicitly the shape of the graphs 
Cfe(s). Thanks to this more visual representation it is 
easily verified that all these graphs represent Od cavities 
of the lattice, because we can place at most one particle 
in any of them. This is a general feature of all FM func- 
tional described by the pattern , so from now on, the 
Cfe(s) will be referred to as Od cavities. 

As we will make clear immediately, the form with 
the Cfe(s) given by Od cavities is a direct consequence of 
the exact dimensional crossover to any Od cavity that FM 
functionals possess. The latter means that if we take a 
density profile which vanishes outside a given Od cavity 
(henceforth a Od profile) and evaluate the functional, we 
will obtain the exact value of the free energy. The only 
known approximate density functionals having this prop- 
erty are FM ones J3, la, |a, llll ■ (As a matter of fact, the 
property can be regarded as the very constructive prin- 
ciple of FM theory Q.) 

Before we start let us define a maximal cavity to be 
any Od cavity which, enlarged by any lattice site, stops 
being a Od cavity because it can accommodate more than 
one particle. Clearly, any Od cavity must be contained in 
a maximal cavity, so Od dimensional crossover needs only 
be proved for maximal cavities. (Notice that there can 
be more than one maximal cavity in a given system.) 

Let us now try to construct the simplest possible func- 
tional of the class which fulfills the exact dimensional 
crossover requirement. Its construction will proceed it- 
eratively. In the first place, if the functional must return 
the exact free energy when evaluated at any Od profile, 
there must appear a term in for each maximal cav- 
ity of the model, and the corresponding coefficient aj- 
must be 1. For the running example we are consider- 
ing, Eq. ||2Jl, this means that we should start off with the 
ansatz 

(3mp] = T.'^o{u), (4) 

because ^ is the only maximal cavity of this model. 



If I0J were the final functional, evaluated at any maxi- 
mal Od profile (one corresponding to a maximal cavity) it 
should return the exact free energy. In the example, all 
maximal Od profiles have the form illustrated in Fig. 
Let us now substitute this profile in and see what 
comes out. For an easy way to do the evaluation, just 
imagine the graphs Cfe(s) as windows which only allow to 
see the content of the lattice nodes they overlap. Then 
the sum over the lattice nodes implies that we must place 
these windows at every lattice site, evaluate the content 
and add up the results of these evaluations. When the 
density profile is a Od one (as in Fig.^), all contributions 
will vanish except those for which the window overlaps 
at least one node of the Od profile. In our example, this 
means that Q will return 

(3Tr [pod] = ^0 it:') + (i-.) + $0 (,^4)] 

i=l 

where pod{s) denotes a Od profile of the form given in 
Fig. (Filled numbered circles in the diagrams rep- 
resent actual evalutations of the density profile for the 
corresponding numbered nodes of the lattice.) 

We can see that, apart from the exact value (the first 
term on the r.h.s.), there appear a number of spurious 
contributions. Therefore Q cannot be the final func- 
tional. These spurious terms are like evaluations of pod 
with non-maximal cavities, so we will try to eliminate 
them by adding new terms to l@J corresponding to non- 
maximal cavities, with the appropriate coefficients. Since 
a term like X]s*^o('>-°) evaluated at podi^) will return 
*o (i»-»2) + ^0 (3H4) + Z^Li (•.), it seems reasonable 
to choose it to remove the first bracket on the r.h.s. of 
(0) (and the vertical one to remove the second bracket). 
Thus we use as our second ansatz 

pmp] = E (o) - *o (-) - $0 m ■ (6) 

Note that in doing this, we have chosen new graphs Cfc(s) 
and their corresponding coefficients ak in When we 
insert pod{s) in this new functional we obtain pj'f^lpod] — 
$0 (^ll^ ~ ^2^=1 ^0 (•;)• We have indeed removed many 
spurious contributions, but there still remain some. It 
should now be clear that in order to remove the latter we 
must add to the previous ansatz the term X^ssz^ ^o(o)- 
This way we obtain the functional (O , which was already 
derived in |3| by a different procedure. It is straightfor- 
ward to check that PT°'^[pod] = ^0 (J^^; thus proving 
its exact dimensional crossover. 

In order to be illustrative, let us apply this procedure 
again to obtain the FM functional for a different model: 
the nearest-neighbor exclusion lattice gas in the trian- 
gular lattice (hard hexagons). This example is different 
from the previous one in that it has two maximal cavi- 
ties: ^ and . Then, the first-step functional must be 
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(7) 



Corresponding to the existence of two different maximal 
cavities there are two different density profiles, as illus- 
trated in Fig. ^p. The exact Od dimensional crossover 
must be satisfied for both of them. Let us start by the 
one with a trianglc-up shape and let us denote it pQd{s)- 
Substituting it in each of the two terms of we obtain 
(using again the window metaphor) 



= *o(,/^,) + 2E*o(..), 

i=l 

= $o(/-^) + $o(%\) + <i>o(i«2) (8) 

+E*o(--)- 



As in the previous example, to remove the "largest" spu- 
rious contributions (those of the dimers) we propose 

nFTip] = /3^r[p] - E + + *o 

(9) 

Substituting pod(s) again we get /3J^r[pod] = *o {j^) - 
Y^i=i ^0 (•>), so we have to add a last correction for the 
point- like cavities, what finally leads to 



/^•^'''[P] = E ['^'0 ( A) + $0 ( V) - *o (/ ) 
-$o('\,)-$o(c«) + $o(°)]- 



(10) 



This functional is exact for pod(s). We would now have to 
check if the same occurs for the Od cavity corresponding 
to the triangle-down in Fig. but symmetry consider- 
ations immediately show that this is the case. In gen- 
eral, checking dimensional crossover for a new Od cavity 
may lead to the appearance of additional spurious con- 
tributions. These have to be eliminated by adding the 
corresponding terms to the functional. 

Finally, notice that H10(l coincides with the functional 
obtained in 0| through a completely different (and far 
more involved) route. 

Let us summarize the procedure to follow for an ar- 
bitrary lattice gas with hard-core interaction. The steps 
are: 

(i) Determine the complete set of maximal Od cavities 
of the model. If we denote them by (fc — 1, . . . , to), 
then the first-step approximation to the functional will 
be [nW(s) = Etec*P(t)] 



(11) 



(ii) Select a maximal cavity and let pod(s) denote a 
generic density profile for it. 

(iii) Insert pod in the current functional, (3J-'['^[p\^ and 
see which spurious contributions appear. Identify the 
terms with the "largest" graphs (those not contained in 
any other of the graphs appearing, except C^) and pick 
one of them. 

(iv) Construct the next step functional l3J'i^i[p\ by 
adding to pj-f^lp] a new term, with its corresponding 
coefficient a^., so that it eliminates the selected spurious 
contribution. 

(v) Repeat steps (iii)-(v) until no spurious contribu- 
tion remains. (Of course, one can exploit the symmetries 
of the model to resume several steps of this process in 
just one, as we have done in the examples.) 

(vi) Repeat steps (ii)-(vi) until exhausting all maximal 
cavities. 

The functional resulting from this process will be of 
the form and will have, by construction, an exact 
Od dimensional crossover. It can be proven that start- 
ing from Hll() there is a unique functional of the form 
with an exact Od dimensional crossover, so any other pro- 
cedure leading to it is equally valid. (In other words, the 
fact that we have chosen to remove the spurious terms in 
decreasing order of "size" is immaterial, but in doing so 
we abbreviate the process.) A sketch of the existence and 
uniqueness proof goes as follows (a more detailed account 
will be reported in 0|). 

Let us form the set V with the lattice C and all max- 
imal cavities C^(s) (s G A; = 1,...,to), and let us 
complete it with all nonempty intersection of any num- 
ber of maximal cavities. For any x,y ^ V, we will say 
that X < y iS all nodes of x are in y. This transforms 

V into a partially ordered set or poset. Any interval 
[x,y] = {z & V : a; < z < y} is a finite subset of V, 
so P is a locally finite poset. Locally finite posets have 
the property [13 that for any mapping f : V V, with 

V a vector space, there exists g : P 1-^ V such that 

fix) = E 9iy), 9ix) = E (12) 

The way to prove this is by inserting the second expres- 
sion into the first, what leads to 



fi{x,x) = l, p.{y,x) 



/i(z,x). 



(13) 



a recursion which defines the (integer) coefficients p{y, x). 
This scheme is referred to in the literature as a Mobius 
inversion, and p,{y,x) is a Mobius function [l7| |. 

For the poset V defined above, let V be the space of 
density functional and take f{x) = the (exact) 

excess free energy functional of a given model on the 
graph X. Specializing H12(l to x — 



(14) 



s£C fc=l 
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where ^'^[p] is an unknown functional. The sum on the 
r.h.s. of H14|l only contains evaluations of J-'^[p\ for Od 
cavities and so is an expression similar to (Q). Now let 
Px'^{s) a generic Od density profile for cavity x. Then, 



ii xDy 

^ynxip] otherwise. 



(15) 



As X n C = X, evaluating H14|l for p^'^{s) yields 



the Bethe approximation [l9j, which is exact in one di- 
mension. Finally, there is a third extension for lattice 
gases in the presence of a porous matrix that we have 
already began to explore [l^- Work along these lines is 
in progress. 

We acknowledge A. Sanchez, C. Rascon and Y. Martf- 
nez-Raton for their valuable suggestions. This work is 
supported by project BFM2003-0180 from Ministerio de 
Ciencia y Tecnologia (Spain). 
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and it is a consequence of Weisner's theorem that * Electronic address: llafuentQmath.ucSm.esI 



vlzjx) = for any x < therefore '^c[Px ] = oi") 
other words, the sum on the r.h.s. of H14(l is exact for any 
Od cavity. 

This completes the proof that the requirement of an 
exact dimensional reduction to Od cavities leads to a func- 
tional of the form Q). As to the uniqueness, it suffices 
to realize that ^{z^x) = (a necessary condition for a 
functional of the form to have an exact dimensional 
reduction to Od cavities) is a particular case of the recur- 
rence (fT^ . whose only solution is p(x,C). 

With this method one can easil y re cover all functionals 
previously obtained in Refs. |J,|a,ll3| and obtain those of 
virtually any other hard-core lattice gas 0] . One further 
striking feature of all functionals obtained in this way 
is that they also have an exact dimensional crossover to 
one dimension, simply because the exact one-dimensional 
functional is of the form Q 0, 0| . 

Clearly the procedure presented above has no restric- 
tion in its application other than the determination of the 
maximal cavities. It can be applied to particles of any 
shape, in any lattice (including regular lattices, Bethe 
lattices, Husimi trees, etc.) and in any dimension. It 
can even be applied to mixtures, either additive or non- 
additive, provided a Od cavity is properly defined as a 
superposition of cavities, one for each species, such that 
at most one particle of only one species can be placed in 
it (see for more details). 

The readers familiar with Kikuchi's cluster variation 
method may have recognized a similarity with the pro- 
cedure we have presented here. The connection is more 
prominent through the Mobius inversion formula 0] and 
will be properly discussed elsewhere 0| . 

Finally, the theory can be generalized in several ways. 
First of all, we have already mentioned that there is a 
straightforward extension to mixtures which recovers the 
functionals for mixtures already derived in 0, 0, Il3| . 
A second extension is the inclusion of "extended" Od 
cavities in which there can be up to n particles. We 
have already checked that the inclusion of two-particle 
Od cavities for the Ising lattice gas (which has repulsive 
and attractive interactions!) yields the functional ob- 
tained from the cluster variation method at the level of 



Electronic address: cuesta@math.uc3m.esl 
[1] M. Nieswand, W. Dieterich and A. Majhofer, Phys. Rev. 
E 47, 718 (1993); M. Nieswand, A. Majhofer and W. 
Dieterich, Phys. Rev. E 48, 2521 (1993); D. Reinel, W. 
Dieterich and A. Majhofer, Phys. Rev. E 50, 4744 (1994). 
[2] S. Prestipino and P. V. Giaquinta, J. Phys.: Condens. 
Matter 15, 3931 (2003); S. Prestipino, J. Phys.: Condens. 
Matter 15, 8065 (2003). 
[3] J. Reinhard, W. Dieterich, P. Maass and H. L. Frisch, 

Phys. Rev. E 61, 422 (2000). 
[4] L. Lafuente and J. A. Cuesta, J. Phys.: Condens. Matter 

14, 12079 (2002). 
[5] J. A. Cuesta and Y. Martinez-Raton, Phys. Rev. Lett. 

78, 3681 (1997). 
[6] L. Lafuente and J. A. Cuesta, Phys. Rev. E 68, 066120 
(2003). 

[7] G. Biroh and M. Mezard, Phys. Rev. Lett. 88, 025501 

(2002) . 

[8] E. Kierlik, P. A. Monson, M. L. Rosinberg, L. Sarkisov 

and G. Tarjus, Phys. Rev. Lett. 87, 055701 (2001). 
[9] S. T. Chui and J. D. Weeks, Phys. Rev. B 23, 2438 

(1981) . 

[10] M. Ya. Azbel, Phys. Rev. A 20, 1671 (1979). 

[11] Y. Rosenfeld, M. Schmidt, H. Lowen and P. Tarazona, 

J. Phys.: Condens. Matter 8, L577 (1996); P. Tarazona 

and Y. Rosenfeld, Phys. Rev. E 55, R4873 (1997); P. 

Tarazona, Phys. Rev. Lett. 84, 694 (2000). 
[12] L. Lafuente and J. A. Cuesta, Phys. Rev. Lett. 89, 

145701 (2003). 

[13] L. Lafuente and J. A. Cuesta, J. Chem. Phys. 119, 10832 

(2003) . 

[14] L. Lafuente and J. A. Cuesta, to be published (2004). 

[15] Applying this procedure, it should now be straightfor- 
ward to obtain, e.g., the FM functional for the first and 
second neighbor exclusion lattice gas in the triangular 
lattice (the so-called t-model of ref. |2) as /9.Fex[p] ~ 

Essz2 [$o (^) + $0 (^) + $0 - 2$o (A) - 

2"l'o ( V) + "^0 (°) ] (^'^^ meaning of the diagrams should 
be self-evident). 

[16] T. Morita, J. Stat. Phys. 59, 819 (1990); Prog. Theor. 
Phys. 115, 27 (1994). 

[17] R. P. Stanley, Enumerative combinatorics, vol. 1 (Cam- 
bridge University Press, Cambridge, 1999), pp. 116-124. 

[18] M. Schmidt, L. Lafuente and J. A. Cuesta, J. Phys.: Con- 
dens. Matter 15, 4695 (2003). 

[19] D. R. Bowman and K. Levin, Phys. Rev. B 25, 3438 

(1982) . 



